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摘要 : 随 着 GPS、GLONASS、Galileo 和 BDS 等 全 球 卫 星 定位 系统 的 建设 和 完善 ， 同 时 
采用 两 个 或 两 个 以 上 定位 系统 信号 进行 多 系统 组 合 定位 成 为 导航 领域 研究 的 热点 。 相 比 传统 
的 单 系统 定位 ， 多 系统 组 合 具 有 更 好 的 定位 、 测 速 、 授 时 性 能 以 及 更 强 的 环境 适应 性 。 通 过 
CAPS、BDS 和 GPS 三 个 系统 的 静态 单 点 联合 定位 解 算 ， 实 现 了 转发 式 定位 系统 与 直 发 式 定 
位 系统 的 联合 定位 ， 验 证 了 中 国 区 域 定位 系统 参与 多 系统 GNSS 定 位 的 可 行 性 ， 并 分 析 了 组 
合 定 位 的 性 能 与 定位 精度 。 分析 表明 , 中国 区 域 定位 系统 参与 多 系统 联合 定位 解 算 是 完全 可 
行 的 ， 单 点 定位 精度 优 于 3m， 且 还 有 提高 的 空间 。 中 国 区 域 定位 系统 的 加 入 ， 丰 富 了 GNSS 
系统 的 选择 , 对 提高 GNSS 联 合 定位 的 多 样 性 、 改 善 联合 定位 性 能 具有 研究 意义 及 应 用 价值 。 
关键 词 : 全 球 导 航 卫 星系 统 ; 中 国 区 域 定位 系统 ; 组 合 定位 
中 图 分 类 号 : P228.1 文 献 标 识 码 :A ”文章 编号 :1672-7673(2018) 


全 球 导 航 卫 星系 统 (Global Navigation Satellite System, GNSS) 经 过 几 十 年 的 飞速 发 展 ， 
美国 的 全 球 定位 系统 (oooaPostonnosysem, GPS)、 俄 罗斯 的 GLONASS、 中 国 的 ssnsspmiasi C Beou navigation 
sateite system, Bos) ~ 欧洲 的 mis Galileo) 均 已 建成 并 投入 使 用 。 在 亚太 区 域内 ， 全 球 定位 系统 单 
点 定位 精度 大 约 为 6-10m， 北 斗 卫 星 导航 系统 由 于 可 见 卫星 数量 较 多 且 卫 星 高 度 角 更 高 ， 单 
点 定位 结果 了 略 优 于 全 球 定位 系统 ， 精 度 大 约 为 4-8ml1。 

为 了 能 够 得 到 比 单 系统 单 点 定位 更 高 的 定位 精度 , 采用 多 个 导航 系统 的 信号 实现 联合 定 

位 导航 是 现今 导航 接收 机 的 一 种 发 展 趋 势 。 相 比 单一 导航 系统 的 定位 , 多 系统 联合 定位 可 用 
的 卫星 数量 显著 增多 , 可 以 通过 星座 优化 选择 几何 布局 更 佳 、 信 号 强度 更 好 的 卫星 实现 组 合 
定位 解 算 ， 能 有 效 提高 定位 的 连续 性 、 准 确 性 和 稳定 性 。 
中 国 区 域 定 位 系统 (Chinese Area Positioning System, CAPS) 是 中 国 科学 院 发 明 的 转发 
式 卫星 导航 定位 系统 ， 是 中 国 二 代 卫 星 导 航 系统 的 重要 组 成 部 分 。 系 统 利用 现 有 的 C 频 段 通 
信 卫 星 组 建 导 航 星 座 , 并 建立 了 一 套 与 导航 相关 的 地 面 设施 。 由 地 面 导 航 主 控 站 的 原子 钟 提 
供 精确 的 时 间 基 准 ， 并 生成 导航 电文 和 测 距 码 ， 上 行 至 通信 卫星 。 通 信 卫 星 上 的 信号 转发 器 
转发 导航 主 控 站 生成 的 导航 信号 ,向 系统 覆盖 区 域内 的 用 户 接收 机 播发 导航 信号 ,使 接收 机 
能 实时 测量 、 解 算 ， 从 而 获得 用 户 接收 终端 的 位 置 、 速 度 和 时 间 信 息 。 为 卫星 定位 方法 的 研 
究 提 供 了 新 的 思路 和 实验 平台 争 。 

本 文 开展 了 转发 式 卫 星 导 航 试验 系统 CAPS 与 BDS、GPS 的 组 合 定 位 解 算 ， 采 用 分 立 的 
BDS/GPS 和 转发 式 CAPS 导 航 接 收 机 分 别 采集 各 系统 的 导航 电文 和 码 伪 距 观测 量 。 通 过 对 三 
系统 卫星 的 优选 构成 导航 星座 , 并 消除 系统 间 钟 差 , 然后 进行 组 合 定位 解 算 。 最 后 通过 统计 ， 
获得 了 有 CAPS 信 号 参与 定位 情况 下 的 组 合 定 位 精度 ， 并 与 BDS/GPS 双 系统 接收 机 定位 结 
进行 了 分 析 比 较 。 
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1 联合 定位 方法 


BDS、GPS 和 CAPS 三 系统 联合 定位 原理 图 如 图 1， 三 系统 的 导航 卫星 共同 组 成 空间 
GNSS 星座 。 由 专用 的 GNSS 接收 机 平台 接收 三 个 系统 全 部 可 见 卫 星 的 导航 电文 ， 测量 码 伪 
距 ， 并 将 原始 定位 数据 保存 到 计算 机 。 最 后 在 计算 机 中 进行 离线 数据 处 理 , 采用 最 小 二 乘法 
进行 位 置 解 算 。 得 到 联合 定位 结果 ， 实 现 三 系统 联合 卫星 定位 ， 并 对 定位 精度 进行 分 析 。 


GNSS 导 航 接收 机 


CAPS 导 航 主 控 站 


图 1 三 系统 卫星 定位 原理 图 


Fig.1 Positioning schematic of three systems 

在 联合 定位 解 算 时 ,CAPS 系统 是 对 传统 卫星 定位 系统 GPS 和 BDS 的 补充 和 增强 .CAPS 
系统 由 于 其 转发 式 定位 的 特色 , 使 用 通信 卫星 实现 定位 。 因 此 可 以 为 多 系统 联合 定位 提供 卫 
星 数量 、 卫 星空 间 布 局 、 灵 活性 和 抗 干扰 等 方面 的 优势 。CAPS 系统 的 加 入 ， 可 以 增加 卫星 
星座 的 有 效 卫 星 数量 ， 优 化 星座 的 空间 布局 ， 减 小 位 置 精度 因子 值 。 此 外 ，CAPS 系统 导航 
信和 号 采用 码 速率 为 10.23Mcps 的 测 距 码 ， 相 比 于 GPS 粗 码 1.023Mcps 和 BDS 的 2.046Mcps 
码 速率 ，CAPS 更 高 的 码 速 率 可 以 得 到 更 精确 的 伪 距 测量 值 。 而且 CAPS 的 系统 时 间 统 一 由 
地 面 站 高 精度 原子 钟 提供 ， 没 有 卫星 间 的 钟 差 ， 也 能 提高 定位 精度 。 


1.1 观测 方程 及 线性 化 


令 卫 星 k 和 接收 机 i 之 间 的 几何 距离 为 p*，c 表 示 光 速 ，dt* 表 示 卫 星 k 的 时 间 相 对 其 系统 
时 之 间 的 时 钟 偏差 ，dt; 表 示 接 收 机 时 钟 偏差 ，e* 表 示 伪 距 的 测量 误差 ,包括 电 离 层 延迟 、 
对 流 层 延迟 、 测 量 噪声 等 , dts 表 示 系 统 间 钟 差 , 包括 dtp 或 dtc, 分 别 为 系统 时 BDT 和 CAPST 
相对 于 GPST 的 时 间 差 。 则 伪 距 友 的 基本 观测 方程 为 : 
L% = pF + c(dtj — dts — dt¥) +e¥, (1) 
其 中 ， 卫 星 和 接收 机 之 间 的 几何 距离 pk 可 由 下 式 计 算 : 


lir 
ut 


EX)? + (Y-Y)? + (Z-Z, (2) 


代入 基本 观测 方程 (0D) 中 ， 得 到 


Lt = J(Xk—X? + (Y*-Y)? + (ZK—Z;)? + c(dt; — dt; — dt") - ef. (3) 
其 中 ， 了 卫星 位 置 (X*,Y*,2Z*) 及 卫星 时 钟 偏差 qt* 可 由 卫星 星 历 获 得 B51]。 
解 算 时 ， 对 于 BDS 和 CAPS 系统 ， 将 接收 机 钟 差 与 系统 钟 差 合并 解 算 ， 即 共有 6 个 未 
知 数 Xi, Yi, Zi dti (dt; 一 dtg)fü(dt; 一 dtc)。 当 组 合 系统 接收 到 信号 的 卫星 数 不 小 于 系统 数 加 
3 时 ( 即 单 系统 至 少 4 颗 卫 星 ， 双 系统 至 少 5 颗 卫 星 ， 三 系统 至 少 6 颗 卫 星 )， 就 可 实现 组 
合 定位 中。 
(3) 式 对 于 接收 机 位 置 (Xi, Yi,2i) 是 非 线 性 的 ， 所 以 在 使 用 最 小 二 乘法 之 前 ， 此 方程 必须 
线性 化 。 式 中 非 线 性 部 分 如 下 : 


f (Xi Y, Zi) = V QC—-X)? + (Y*-Y)? + Z-Z)^, (4) 


需要 选择 一 个 接收 机 的 初始 位 置 (Xio, Yo, Zao) 开 始 线性 化 ， 此 处 选择 初始 位 置 为 地 心 
(0,0,0). XE SOS EAX, AY, AZ 如 下 ， 用 于 更 新 近似 的 接收 机 坐标 ，?Z 为 迭代 计算 次 数 ， 每 次 旭 
代 计 算 均 以 前 一 次 计算 更 新 后 的 位 置 作为 本 次 计算 的 初始 位 置 : 

Xip+1 = Xip + AXi 
Yip+1 = Yip t AY; p = 0,1,2…, (5) 
Zip+1 = Zip + AZi 

将 (4) 式 在 每 次 迭代 计算 的 初始 位 置 (Xi,, Yip Zip ) EAS ETT ,— 8: 

Of CY Tip) iy 十 Of (Xip Yip Zip) ay 


0X Lp ip 
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4 Of Mp Yip Zip) (tig Yip i2) hy, (6) 


fis Yip+1 Zip+1) m f (Xiy, Yip Zip) + 
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Of (Xi Yi, Zip) E (XF-Xiy 


OXip pr 
Of (Xip Yip Zip) = Y*-Y,, (7) 
OYip pi 
Of (Xip Yu, Zip) " ZZ 
0Zip u m 
至 此 ，(3) 式 变 为 一 阶 线性 观测 方程 ; 
k k k 
VE = pl, - P AX, — — L2 AY, — L p AZ, + c(dt; — dts — dt*) + ek (8) 
Pip Pip Pip 


其 中 ，p 包 为 由 接收 机 近似 位 置 解 算 的 距离 : 


(9) 


1.2 最 小 二 乘法 


最 小 二 乘法 通常 用 于 解决 线性 方程 组 中 方程 数 大 于 未 知 数 时 ， 方程 无 解 的 情况 Bi。 对 于 
一 个 无 解 的 方程 4x = b， 其 中 4 有 m 行 nh 列 ，m > n。 即 方程 中 观测 值 bj,… ,bx 的 个 数 大 于 参 
数 x1,…, xn 的 个 数 。 此 方程 误差 向 量 为 6 = b — AR, 其 平方 值 ||lell? = (b — Ax) T (b 一 Ax) 是 m 
个 分 立 误差 的 平方 和 ， 为 得 到 使 误差 向 量 长 度 最 小 的 最 优 解 介 ， 最 小 化 该 二 次 方程 得 到 : 


ATA$ = ATh 或 ® = (ATA) 1ATb . (10) 
于 是 ，(8) 式 可 重新 写 为 矢量 形式 ， 即 


AX; 
XE—X, Y*—y, Z*-Z AY; 
kpk p ip ip _ ^ ip i ode. de 
Li = Pip | p, PE, PE, 1 AZ; cdt* t ej. 
í c(dt, — dt;) 
(11) 


再 将 其 整理 成 最 小 二 乘 问题 4x = b 的 常规 表达 式 ， 即 
X*-X;;  Y*-Yi 2*-Zip 
| Pip Pip Pip 
从 以 上 单一 方程 中 得 不 到 唯一 解 ， 故 令 bk = LE— pi, 十 cdt* 一 ek， 联 立 观测 的 各 导航 系统 
卫星 方程 组 成 方程 组 ， 则 最 终 解 可 由 (13) 式 获得 ， 其 中 系数 矩阵 4 的 前 mm 行 对 应 为 GPS 卫星 ， 
(X56,Y6,26) 为 GPS 卫 星 坐 标 ， 中 间 n 行 对 应 为 BDS 卫 星 ，(X3,Y3,23) 为 BDS 卫 星 坐 标 ; 最 后 
! 行 对 应 CAPS 卫 星 ，(X5,YC,2Z5) 为 CAPS 卫 星 坐 标 。 


1 Xy, = LE — pi, + cdt* — ef. (12) 
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"m -nclz68[, HUE—fE: AXip+1AYip+1AZiw+1。 此 解 可 以 车 加 到 接收 机 的 初始 
近似 位 置 上 ， 从 而 得 到 下 一 步 更 精确 的 近似 位 置 : 
Xip+1 = Xip t AXip+1 
Yip+1 x Yip + AY; p+1 . (14) 
Zip+1 = Zip + AZip+1 
H AEH X pro Yipep Zip INE Xip Yip Zip» 代入 (13) 式 循环 计算 ， 直 至 q 次 之 后 方 
程 的 解 AXiy4o,AYi s+, AZip+a 达 到 期 望 的 误差 数量 级 ，q 即 为 最 小 二 乘 先 代 次 数 。 通 常 ， 四 
到 五 次 循环 计算 后 便 可 符合 目标 。 


2 ”试验 方案 及 过 程 


2.1 试验 方案 


基于 BDS, GPS 和 CAPS 的 组 合 卫星 定位 试验 采用 静态 单 点 定位 方式 。 试 验 接收 机 安 
装 在 中 国 科 学 院 国家 天 文 台 楼 项 。 由 于 CAPS 系统 工作 在 C 频段 ， 与 其 他 两 系统 采用 的 L 
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频段 不 同 ,为 了 改善 信号 接收 性 能 ， 采 用 分 立 的 接收 天 线 靠近 放置 。 导 航 射 频 信号 转换 到 中 
频 后 ， 分 别 由 CAPS 接收 机 和 BDS/GPS 双 系 统 单 频 接收 机 基带 捕获 跟踪 解 调 ， 并 实时 同步 
记录 三 系统 的 卫星 星 历 、 伪 踊 等 原始 定位 数据 ， 并 存储 到 计算 机 。 最 后 在 计算 机 中 采用 
MATLAB 软件 编程 进行 数据 处 理 和 联合 定位 解 算 ， 并 分 析 试 验 结果 ， 统 计 定 位 精度 。 定 位 


试验 系统 框图 如 图 2。 
CAPS GPS/BDS 
B SERE 
LNA LNA 


y Y 
C-IF < 一 一 频 综 一 > LIF 


i | i 


导航 基带 1 < 一 一 原子 钟 — 导航 基带 2 


试验 工控 机 
原始 定位 数据 1 联合 定位 解 算 原始 定位 数据 2 
数据 分 析 


4] 2 三 系统 卫星 定位 试验 原理 框图 


Fig.2 Schematic of combined positioning experiment 


2.2 试验 过 程 


(1) 在 天 文 台 楼 项 放置 一 台 BDS/GPS 双 系 统 定位 接收 机 ， 进 行 48 h 连续 测量 ， 取 其 定 
位 解 算 结果 数值 作 平 均 ， 获 得 静态 定位 参考 点 坐标 值 (xo, yo, Zo) 73(-2173764.5307, 
4383137.0366, 4078283.8735) 

se A er i a ea rd 
包括 全 部 可 见 GPS, BDS 和 CAPS 卫星 的 卫星 星 历 、 卫 星 钟 差 、 电 文 修正 信息 、 码 
伪 距 测量 值 等 。 

(3) 选取 各 系统 同时 段 的 卫星 定位 信息 ， 由 卫星 星 历 计算 每 个 历 元 的 卫星 位 置 坐标 值 ， 
同时 也 采用 三 系统 全 部 可 观测 到 的 卫星 用 最 小 二 乘法 进行 定位 解 算 。 

(4) 处 理 定位 结果 数据 ， 分 析 统 计 联 合 定位 时 定位 结果 的 准确 度 和 精度 ， 与 BDS/GPS 
双 系 统 定位 结果 作对 比 ， 得 出 试验 结论 。 


2. 3 数据 统计 方法 


当 设 备 处 于 正常 工作 状态 ， 三 系统 均 成 功 捕获 跟踪 信号 后 ， 连 续 记 录 m (n > 1000) 个 历 
元 定位 伪 距 测量 数据 及 该 期 间 的 导航 电文 。 对 此 n 个 历 元 数据 进行 定位 解 算 ， 得 到 每 历 元 的 
定位 解 算 结果 。 采 用 的 计算 步骤 如 下 : 

(1) 对 每 个 历 元 ， 计 算 名 个 方向 的 定位 偏差 及 三 维 合成 总 误差 : 


Az; = Zi — Zo i212,-,n (15) 


Api = |Ax;? 十 Ay;? 十 Az? 
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其 中 : 
Xp Yp Zi 为 第 i 个 历 元 定位 结果 的 三 维 华 标 ， 单 位 为 ms 
Xo Yo» Z0 为 静态 定位 参考 点 坐标 ， 单 位 为 mi 
(2) 分 别 计算 各 个 方向 及 三 维 合成 总 误差 的 误差 平均 值 和 误差 标准 差 : 


n 
—! 
pl) (16) 


1 
n—1 


X n-m ， (17) 
i=1 


ZN 中 : 
Ap 为 n 个 历 元 三 维 合成 总 误差 的 误差 平均 值 ， 单 位 为 m; 
Sap 为 n 个 历 元 三 维 合成 总 误差 的 误差 标准 差 ， 单位 为 m; 


3 ”试验 结果 


试验 选取 2017 年 9 月 11 日 18 时 19 分 至 9 月 12 日 10 时 34 分 的 测量 数据 ;使 用 了 全 
部 可 见 的 GPS、BDS 和 CAPS 卫星 数据 , 每 一 时 刻 能 同时 接收 到 的 卫星 最 少 为 16 9i, CAPS 
系统 所 用 卫星 详情 如 表 1。 
表 1 CAPS 系统 用 于 本 试验 的 卫星 情况 表 
Table 1Satellites of CAPS used in experiment 


卫星 编号 卫星 类 型 卫星 轨 位 
5 SIGSO 148°E 
6 GEO 110.5°E 
8 GEO 76.5"E 


采用 MATLAB 编制 CAPS/BDS/CAPS 多 系统 解 算 程 序 进 行 位 置 解 算 ， 解 算 结 果 画 出 平 
面 散 点 图 如 图 3。 
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图 3 双 系 统 / 三 系统 定位 解 算 误差 平面 散 点 图 
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Fig.3Positioning horizontal error of GPS/BDS and CAPS/GPS/BDS 


双 系 统 、 三 系统 组 合 时 ， 


位 置 精 度 因子 值 情 况 如 图 4。 
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图 4 双 系 统 / 三 系统 定位 解 算 卫星 PDOP 值 变化 图 
Fig.4PDOPvalues of GPS/BDS and CAPS/GPS/BDS 
定位 精度 统计 结果 见 表 2. 
表 2 双 系 统 /三 系统 系统 定位 精度 统计 表 
Table 2Positioning accuracy of 2/3 systems 
定位 系统 BDS+GPS BDS+GPS+CAPS 
水 平定 位 精度 /m 1.2534 1.1270 
高 程 定位 精度 /m 1.6842 1.7603 
三 维 定位 精度 /m 2.0994 2.0902 


4 ”结论 


态 单 点 定位 情况 下 ，CAPS 系统 的 加 入 ,提供 了 更 丰富 的 卫 


星 选择 ， 改 善 了 空间 卫星 星座 构成 ， 减 小 了 位 置 精 度 因 子 值 。 因 此 提高 了 定位 精度 ， 尤 其 是 


水 平定 位 精度 提高 较为 明显 。 在 以 后 的 研究 中 可 以 通过 进一步 精细 消除 误差 或 采用 载波 相位 


[三 


测量 值 参与 解 算 ， 定 位 精度 还 


有 较 大 的 提高 空间 。 


Analysis of combined positioning accuracy based on CAPS, 


BDS and 


GPS 


Fu Shengyou"^, LiShengming’, Wang Zhaoruji 


(1. University of Chinese Academy of Sciences, School of Astronomy and Space Science, Beijing100049, China) 


(2. Chinese Academy of Sciences, National Astronomical Observatories of China, Beijing100012, China) 


Abstract: Withconstructing and completingof global satellite positioning systems such as 


GPS, GLONASS, Galileo and BDS, using two or more than two signals ofdifferent positioning 


systems toachieve the combined positioning has become a hotspot in research of 
navigation.Compared with traditional single system positioning, multi-system combination has 
stronger environmental adaptability and better performance of determining position, velocity and 
time.By combined static single point positioning of CAPS, BDS and GPS, we achieve the 
combined positioning of transponding system with straight forward system, verify the feasibility 
of CAPS participating in GNSS  multi-system positioning, and analyze the positioning 
performance and accuracy.The analysis shows that CAPS are fully feasible in multi-system 
positioning, single point positioning accuracy is better than 3 meters and there is room for 
improvement.The addition of CAPS enriches the selection of GNSS system, and has research 
significance and application value for improving the diversity of GNSS positioning and the 
positioning performance. 
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